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Abstract 

Variational principles are developed within the framework of a spinor representation of the surface 
geometry to examine the equilibrium properties of a membrane or interface. This is a far-reaching 
generalization of the Weierstrass-Enneper representation for minimal surfaces, introduced by mathe- 
maticians in the nineties, permitting the relaxation of the vanishing mean curvature constraint. In 
this representation the surface geometry is described by a spinor field, satisfying a two-dimensional 
Dirac equation, coupled through a potential associated with the mean curvature. As an application, 
the mesoscopic model for a fluid membrane as a surface described by the Canham-Helfrich energy 
quadratic in the mean curvature is examined. An explicit construction is provided of the conserved 
complex- valued stress tensor characterizing this surface. 

1 Introduction 

The coarse-grained description of a physical surface is often entirely insensitive to microscopic details, 
involving an energy function that depends only on the surface geometry. A familiar example is provided 
by a soap film or an interface which will reach equilibrium when its surface area is minimized. In the 
absence of a pressure difference across the surface, this state is described by a minimal surface characterized 
by a vanishing mean curvature; it will thus assume a saddle shape almost everywhere. 

From a mathematical point of view minimal surfaces are very special: the functions embedding the 
surface in space are harmonic so that, in an appropriate parametrization, the surfaces will be described by 
analytic functions. This is the essential element of the classical Weierstrass-Enneper (WE) representation 
of minimal surfaces, introduced 150 years ago, which has played a central role in the development of the 
subject ever since [TJ|2]. While physically it is often adequate to settle for a small gradient approximation 
in terms of a height function, where one can get by with a more modest toolkit, the WE representation 
of the surface geometry provides an indispensable handle on the global features of minimal surfaces; the 
Schwartz P Surface or plumber's nightmare [3], describing liquid crystalline phases of fluid membranes, 
provides a good case in point [3]. 

Even if one is interested in interfaces, however, one almost always needs to accommodate a pressure 
difference across this surface; the minimal surface is replaced by one of constant mean curvature. Never 
mind addressing questions of stability or examining fluctuations, simple as this modification may be, its 
description lies beyond the scope of the WE representation. 



There is, however, a generalization of the WE framework that applies to any surface. Back in 1979 Ken- 
motsu showed how the classical representation could be tweaked to describe surfaces with any prescribed 
mean curvature [5]; a decade later, the WE representation was reformulated in terms of a two-component 
spinor field, which also turned out to be its natural setting By the mid- nineties, this framework had 

been extended to accommodate a non- vanishing mean curvature [El [S] ■ As emphasized by Konopelchenko 
and Taimanov in |10j . it was no longer even appropriate to think of the mean curvature as prescribed. 
For reviews see references [TT] and |12) . 
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The spinor appearing in this framework satisfies a linear Dirac-type equation on the complex plane; its 
components couple through a real-valued external potential. As in the classical WE representation, the 
functions which describe the embedding of the surface in space occur as integrals of closed differential 
one-forms quadratic in the spinor field. The mean curvature of the surface is proportional to the poten- 
tial. When it vanishes, the components uncouple and the original representation for minimal surfaces is 
recovered. 

Our focus will be on the mathematical description of fluid membranes on mesoscopic scales [13] . On these 
scales the membrane is described with unusual accuracy as a two dimensional surface; the energy associated 
with a given configuration is the Canham-Helfrich bending energy quadratic in the mean curvature [14] . 
Unlike any familiar elastic material, in-plane shear goes unpenalized; thus the surface is characterized 
completely by its geometrical degrees of freedom and membrane elasticity can be addressed in terms 
of geometry. A striking feature of this energy in the spinor representation is that the integrated mean 
curvature squared depends only of the potential. Modulo a topological contribution to the energy, this is 
also the most general geometrical energy associated with bending. It is thus curious that there has been 
no systematic attempt made to look more closely at this functional in terms of these variables. In this 
paper we will examine the behavior of the Canham-Helfrich energy under variations of the potential; in 
particular we will show how the equation describing equilibrium shapes-non-linear both in the spinor and 
in the potential-is obtained. 

What would appear to be the obvious way to do this also turns out to be riven with technical difficulties: 
if the spinors and the potential are treated as the fundamental variables, the relationship connecting their 
variations to that of the embedding functions defining the surface is non-local. One way to sidestep this 
difficulty is to introduce an appropriate set of Lagrange multipliers, a spinor analogue of a framework 
introduced a few years ago to reformulate the variational principles for parametrized surfaces in terms 
of the induced metric and the extrinsic curvature [15) . Just as the metric and the extrinsic curvature 
are not independent variables, the spinor and the potential are also constrained (by the Dirac equation). 
This is just as well; for a naive variation of the potential appears to imply that the only critical points of 
bending energy are minimal surfaces, a manifestly incorrect conclusion. In this framework, the relationship 
connecting the spinors to a surface, characterized by its three embedding functions, will also be input as a 
constraint. The multiplier enforcing this constraint quantifies the local change in energy as the surface is 
deformed; it is thus identified as a "stress tensor" . Unlike the stress in a typical elastic medium, however, 
it will be determined completely by the geometry. A surface in mechanical equilibrium is described by a 
conserved stress tensor. 

The description of equilibrium in terms of a geometrical stress tensor has been discussed previously in the 
context of a parametrized description of the surface; in reference [16j it was identified as the Noether current 
associated with the translational invariance of the energy, an approach refined in reference [15] using the 
framework of auxiliary variables. More physical continuum mechanical/thermodynamical treatments have 
been presented in references [T7| and [THJ ). Work in the latter direction dates back to the work of Evans 
[19j . The concept of geometrical surface stress is also familiar, if only implicitly, in the context of minimal 
surfaces [2Ql [21] where they appear in the identification of the weights-external forces-that provide the 
tension preventing the collapse of the surface. 

A problem of significant current biophysical interest is the study of interactions between surface-bound 
particles mediated by the surface geometry, with the particles in question identified as proteins (see 
reference [22] . for example, and more recently [23). In [24 this problem was approached in the non- linear 
regime by examining the geometrical stresses associated with the deformed geometry. By casting the 
Euler-Lagrange equations as a conservation law it became possible to gain access to global information 
which is not manifest when the divergence is dismantled into tangential and normal parts. Just as the 
techniques of complex variables can be exploited very effectively in the linear regime (as demonstrated 
spectacularly in reference [22] ). the spinor framework renders it possible to apply these same tools to 
examine physical processes which lie beyond the scope of a linear description. Recent reviews, which 
summarize nicely the role of geometrical stresses in a soft matter and biophysical context, have been 
provided by Deserno [35] and Powers [2"S] . 

The simplification of the surface geometry in this representation involves giving up, to a large degree, the 
reparametrization invariance inherent in the geometrical nature of the problem. An arbitrary deformation 
of the surface will not be compatible with the manifestly conformally flat form of the metric. For example, 
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any purely normal deformation of a non-planar surface will generally be inconsistent with the represen- 
tation; consistency will require a compensating tangential deformation to tag along. Fortunately, this 
awkward technical point never needs to be addressed explicitly in the variational framework we introduce. 
It does, however, manifest itself in the conservation law for the stress tensor. Whereas the tangential con- 
servation laws would be expected to be satisfied identically in any completely reparametrization invariant 
framework, where small tangential deformations of the geometrical degrees of freedom are identified with 
reparametrizations (see [27]), here they are not. This is also, of course, just as well: the Euler-Lagrange 
equations for the spinor and the potential leave undetermined the off-diagonal component of the tangential 
stress; the tangential components of the conservation law turn out to provide the differential equations 
determining this missing component. 

It is evident that the description of equilibrium in terms of a conserved stress tensor possesses a number 
of interesting model-independent properties that are best appreciated by extending the discussion to 
accommodate "energy" functionals of a more general kind. We will thus examine geometrical energies 
involving powers of the mean curvature. The Canham-Helfrich model occurs as a special case. 

We begin in section [2] with a brief review of relevant features of the spinor representation of surfaces. 
Our main results are contained in section [3] where variational principles are developed within the spinor 
framework. In particular, the energy most relevant in condensed matter applications, consisting of a linear 
combination of area, the integrated mean curvature and bending energy is considered. We conclude with 
a brief discussion as well as a few suggestions for future work. Various useful identities are collected in a 
set of appendices. 



2 Dirac equation for surface spinors 

In the generalized WE representation of the surface E the three functions describing its embedding in 
three-dimensional Euclidean space E 3 are expressed in terms of a 2-component spinor field ip(z, z) = 
(-01 (z, z), ip2(z, z)) defined on a simply connected domain D of the complex plane C (the overbar denotes 
the complex conjugate). This spinor is a solution of the Dirac equation 

Vip = 0, (1) 

where T> is the first order differential operator, 

»-U*)+(i;s). < 2 > 

involving a real valued potential V. This potential couples the spinor components; Eq.([T]) reads 

d,ih = -vih, dgih = vih. (3) 

The spinor £ defined by £ = (— V>2, V'i) T also satisfies the Dirac equation: 

Vi = . (4) 

We will also write £ = where 

c-O". 1 )- m 

Note also that £ unlike ip transforms like ip under the action of SU{2). A linear combination of tp and £ 
corresponds to a rotation of the surface, as described in | Appendix A| 

The immersion of D into E 3 defined by ip(z,z) is given by 

X(z,z) = / 4>(w,w) , (6) 
J 1 

where (f>(z, z) is the vector- valued 1-form with the following component^] 

^(z,z) = Re[{^ 2 2 -Vi 2 H, (7a) 

2 (z,z) = Re li^+^dz] , (7b) 

(/) 3 (z,z) = 2Re[in<ip2dz] ; (7c) 

1 Here we interchange cf) 1 and <f> 2 with respect to the more usual definitions in the literature, for example 1 1 21 . as this 
facilitates the connection with the classical WE representation. The number of minus signs appearing in the calculations is 
also reduced. 
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7 is a path on the complex plane terminating at the point z. It is evident from their definition that the 
embedding functions are invariant under complex conjugation, ip — > ip. 

It is simple to show that the 1-forms (JTJ) are closed in the complex plane, i.e. dip 1 — 0: one uses the 
fact that the spinor ip satisfies the Dirac equation, ([T]) and that V is real. Note that the derivatives 
d z ipi and dzipz never need to be evaluated in this argument; in addition, it should be remarked that the 
corresponding 1-forms constructed with the imaginary part in Eq.((7]) are not closed. 

Closure and Stokes theorem together imply that the definition of the embedding functions is indepen- 
dent of the choice of the path 7 and thus well defined. It is worth pointing out that the Dirac equation 
(JT|) appears to be the most general linear equation for the spinor consistent with closure. 

2.1 Intrinsic geometry of the surface 

The tangent vectors to the surface adapted to this parametrization are e 2 = 9 2 X and e 2 = e 2 ; in terms 
of the spinor components e 2 is given by 

! / W-^i 2 \ 
e, = ^ i(t> 2 2 +^ 2 ) • (8) 

The vector- valued 1-form </> defined in Eq.© can thus be expressed in the alternative form <f> = 2Re (e z dz) . 

The two tangent vectors are null with respect to the scalar product in E 3 (denoted by •), in other words, 
their norms vanish: e 2 • e 2 = = e 2 ■ e 2 . However they are not orthogonal: e 2 • e 2 = 1/2 \tp\ , where 
|^| 2 = ^ = |Vi| 2 + |^2| 2 and|^| 2 = ^i- 

This parametrization is isothermal (or conformal): the line element is given by 

ds 2 = HVl 2 ; (9) 

the induced metric on the surface assumes the form 

wvo n (10) 



9ab- 2 y t Q 

So the metric tensor is a Weyl rescaling of its Euclidean counterpart on the complex plane, i.e., it is 
manifestly conformally flat. The fourth power of the spinor norm provides the conformal factor. The 
intrinsic geometry of the surface is completely determined by \ip\. 

The residual reparametrization freedom consistent with the isothermal form of the metric is captured by 
analytic functions. Under the holomorphic transformation z — > w(z), z — > w(z), one finds that 

-01 (z, z) -> u/(2:) 1/2 V>i(w, w) , ip 2 (z,z) -> w'(z) 1/2 ip 2 (w,w) , V(z,z) -> \w'(z)\V(w,w) . (11) 

These variables are thus scalar densities. 

The determinant of the metric is g — — 1/4| -0 1 8 , so that ^fg — i/2\ip\ 4 is imaginary. The area is given by 

Despite appearances, it is a real form: dA = —i/2\ip\ 4 dz A dz = dA. 
2.2 Extrinsic geometry of the surface 

The normal vector to the surface E defined by n = e 2 x e s /^/g, in terms of the spinor components is given 
by 

2 Re (V>i^2) 

n=^| 2/m(V^ 2 ) ) . (13) 

\H 2 -\H 2 
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The second fundamental form or extrinsic curvature tensor denned by K a b = — n • d a eb, is given by 

K h -( A V ^ 2 ) (14) 
Kab ~ { V\^\ 2 A ) ' (14) 

where A denotes the Wronskian of ipi and $2, A = ^dz^Pi — A is invariant under rotation, as 

are the norm of ip and V. Under reparametrization, z — > w(z), A transforms as a scalar density: 

A{z,z) -> w'(z) 2 A. (15) 

Besides the potential and the spinor norm, the Wronskian is the other function of the spinor and its 
derivatives which occurs naturally in this framework. 



The eigenvalues of the shape operator K a \> — g ac K c b (g ab is the inverse metric) are the two principal 
curvatures, C\ and Ci- The two symmetric curvature invariants, (twice) the mean curvature K = C\ + C2 
and the Gaussian curvature Kg = C1C2 are given respectively by 

* = ]^V, (16) 
Kg = ^(V 2 W-\A\ 2 ). (17) 

If V = the surface is minimal (K = 0). By defining the components of the spinor ip in terms of two 
analytic functions f{z) and g{z) as ipi = j x l 2 gl\[2 and ip2 = f 1 ^ 2 /V% we obtain the <fi corresponding to 
the original WE representation for minimal surfaces [5] given by 



<p = Re 



f -(l-g 2 ,l + g 2 ,2g) T dz 



(18) 



The classical Gauss-Codazzi integrability condition on the surface geometry ( |Appcndix D I identifies 
(twice) the Gaussian curvature (|17p with the intrinsically defined Ricci scalar given by (IB. 41) . Thus 
\A\ is completely determined once V and are known; it is a measure of the difference of the principal 
curvatures: 



\A\*= {Cx -^ ? W. (19) 
The phase of A captures the principal directions; A vanishes only at umbilical points. 

2.3 Spinor presentation of bending energy 

Bending energy is quadratic in curvature. For a two-dimensional surface this implies that it must also 
be scale invariant. There are several different possibilities. The simplest of these is positive definite, 
proportional to 

Hi = \ JdAK ab K ab = !JdA(C 1 2 + C 2 ) , (20) 
which vanishes on a planar region with C\ — = C2 (K a b = 0). The Canham-Helfrich energy [14] 

H 2 = \ J dAK 2 =4i J dzAdzV 2 , (21) 

which vanishes on a minimal surface, assumes a remarkably simple form in the spinor representation. 
The conformally invariant Willmore energy is given by the integrated squared difference of the principal 
curvatures [55] 

H 3 = [ dAK ab K ab = i / ' dA{d- C 2 ) 2 , (22) 



2 . 

where K ab — K ab — \g a bK is the trace-free curvature. H3 vanishes on a sphere. It is simple to see that 

H 3 =4i [ dzA dz\-prr , (23) 
J W A 

involving the modulus of the Wronskian A. 
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Finally, one has the Gaussian bending energy, H\ = J cLAKq which is identified as the Gauss-Bonnet 
topological invariant. It is clear that any quadratic bending energies can be constructed using any two 
of Hi, H2, H3 and i/4. Choosing one of these as H4, it is clear that any one of the remaining three 
is equivalent to any other modulo a boundary term associated with the Gauss-Bonnet invariant. In 
particular, the Willmore energy differs from the Canham-Helfrich bending energy by a topological term 
linear in the Gaussian curvature. Thus, despite the functional differences between the two, their Euler- 
Lagrange equations had better be identical. As a consistency check, we will confirm this fact explicitly in 
the spinor framework in |Appcndix F| 



3 The construction of the stress tensor 

We are interested in examining the variational properties of various physically relevant energy functionals 
of the surface geometry within the spinor framework. This may be the area, bending energy or some other 
geometrical invariant. Such an invariant will generally be some scalar density, L, constructed using the 
spinor and the real potential V, integrated over the the surfaced 



H = iJ dz A dz L(ipi,ipi, ip2, 1P2, V) . (24) 

Derivatives may also appear in L. Derivatives of ip, if they occur, are expressible in terms of ip itself and 
derivatives of ln|?/>| 2 , V, and A. Where higher derivatives are involved it is not entirely straightforward 
in this framework to identify scalar combinations of the variables appearing in the argument of L; it is 
facilitated, however, by taking into account the transformation properties of the densities ip, V given by 
Eq. (|ll[) and of A given by Eq. (|15p . and forming products with weight zero. 

There are two sets of constraints that need to be accommodated when H is varied with respect to these 
variables: the Dirac equation ([1} implies that the variations in the spinor and the potential are not 
independent; we also must encode how to reconstruct the surface in terms of the spinors. These constraints 
are enforced in the variational principle by introducing a number of Lagrange multipliers. Thus we 
construct the following functional: 

H c = H + i J dz Adz X^Vi/j + i J dzAdz A T X>£ 

+ if dz A dzf z ■ (e z - d z X) +i J dz Adzf z ■ {e s - d z X) , (25) 

where V is defined by Eq.@ and e 2 by Eq.©. The Lagrange multipliers A and P enforce respectively 
the Dirac equation and the identification of the tangent vector^]. It is now legitimate to treat X, t/j, 
■0 and V as independent variables. The multiplier A is a 2-component spinor defined by A = (A 1 , A 2 ) 

@. The multipliers P form a vector in E 3 defined by P = (Z 1 ,/ 2 ,/ 3 ) 7 - This construction is a spinor 
counterpart of that for parameterized surfaces [15] . There are, however, a number of important differences 
in its implementation. 

It may appear, at first, that the tangency constraints are redundant is the spinor framework. After all, 
any spinor satisfying the Dirac equation, will define a surface, and a naive counting of degrees of freedom 
suggests that nothing is amiss if they are dropped. This would certainly simplify the variational principle; 
unfortunately, it would also be wrong. The multipliers P, as we will see, will get identified with the stress 
tensor in the surface. Failure to enforce the constraint would leads one to the invalid conclusion that 
equilibrium solutions always have vanishing stress which is clearly not the case. 

In components, the functional H c assumes the for 

H c = H + i J dzAdzX 1 (d z ip 2 + Vtpi) +i J dzAdzX 2 (d z tpi - Vip 2 ) + c.c. 

+ ijdzAdz(t- (j 2 2 ~ <Ai 2 ) + i Y (W + ^1) + - P ■ d z X^j + c.c. . (26) 

2 It will be convenient to work with a scalar density L, rather than a scalar C; the two are related by L = jIV'l 4 ^-- 
3 As defined, both A and f z are densities. 

4 The fact that that the Dirac equation Q is the complex conjugate of JTJ is reflected in the fact that the multipliers 
enforcing these conditions are also complex conjugates. In addition, f z = f z . 
5 c.c. represents the complex conjugate expression. 
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Begin with the variation of the embedding functions X. Performing integration by parts to collect the 
derivatives of f z and f z in a divergence, one finds 



5yiH c =i J dzAdz (dj z ■ SX - d z (f z -SX)) + c.c. . 

Thus, the Euler-Lagrange derivative with respect to X, ex = 5H C /SX is identified as 

ex = dj z + dd s ■ (27) 

Critical points of H c satisfy the Euler-Lagrange equations, ex = 0. The solutions of the Euler-Lagrange 
equations are therefore described in terms of the conserved complex- valued surface "stress tensor" , f z |£|. 
It remains to construct f z explicitly. This will involve the solution of the Euler-Lagrange equations for 
the variables tp, ip and V. 

Varying He with respect to the potential V gives 

S V H C = ijdzAdz(^ + AVi + A 1 ^! - A 2 V>2 - X 2 ^ SV ; 

the corresponding Euler-Lagrange equation is thus 

e v = AW + X T a^ + ttt = . (28) 
ov 

Here 03 is the Pauli matrix with 1 and — 1 along the diagonal. Since V is real, this equation is also real. 
The spinor A thus satisfies a single inhomogeneous linear algebraic equation in a two-dimensional complex 
vector space. One solution of this equation is given by 

A * = -2W^- (29) 

However, this solution is not unique. Let us write Ea. (|2"5f in the form Re(X^a3tp) = —5L/SV. Any 
spinor A with an imaginary projection onto a^ip is evidently a solution of the homogeneous equation 
Re(\^<j3ip) — 0. This equation has two solutions. The first of these is of the form 

X r h = ira 3 Tp, (30) 

where r is an arbitrary real-valued function. This is because X^a^ip — ~ i\'4'\ 2 i' which is manifestly 
imaginary. The second solution, which involves the spinor £, is given by 

X c h = ca 3 £, (31) 

where c is an arbitrary complex- valued function. This solution satisfies X^a^ = on account of the 
orthogonality of tp and £. 



A = ( -TTTTFlS+H^ +«t 3 £, (32) 



Thus the complete solution A = A; + Xh is given by 

1 5L 
~2\ij\ 2 6V 

the components of which read 

Al = (" 2WS + ^ " C "^' A ' = " V ' 2 " C_ ^ 1 ' (33) 

This solution possesses 3 degrees of freedom per point: one for r and two for c. To justify the counting, note 
that in terms of real- valued variables, the solution of equation (|28)) describes a 3-dimensional hyperplane 
in a 4-dimensional vector space equipped with a (non-degenerate) inner product of signature (+, +, — , — ). 
It will be shown below that the ambiguity reflected in the functions r and c is a gauge artifact associated 
with the parametrization. 



6 With a densitized stress tensor covariant differentiation is replaced by partial differentiation. 
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Let us introduce the quantities /+ = f 1 + if 2 and / = f 1 — if 2 . The Euler-Lagrange equations for the 
spinor ip are given b}|3 

e*! = ~Af~ +M 3 +T 1 =0, (34a) 
= V2/++^i/ 3 +T 2 = 0, (34b) 

where 

T 1 = -^--0 s A 3 + A 1 V, T^-^-fl^-A'V. (35) 
together with their complex conjugate counterparts. 

Solving this set of equations for f 2 is facilitated by by first expressing the Cartesian components / + , /~ 
and / 3 in terms of their geometrically more relevant counterparts with respect to the basis of tangent 
vectors adapted to the surface {e z ,e 2 ,n}. The latter decomposition of f z is given by 

y/gf = i (f 2 ■ e 2 e z + P ■ e z e 2 ) + y/gf z • n n . (36) 

The three projections can be expressed in terms of the Cartesian components, / + , / - , f 3 , as follows 

fz = P ■e z = ±(4, 2 2 f + -^ 2 r)+M2f\ (37a) 

f% = f .e z - = -i(^ 2 / + -^ 2 2 r)+^i^/ 3 , (37b) 

f = f*-n=j^p + ViV*/ - + (l^i| 2 -|^ 2 | 2 )/ 3 ) ; (37c) 
these three equations are now inverted in favor of / + , / _ and / 3 to obtain 

/ + = (V- 2 2 n - ^i 2 A + l^lVi^/*) , (38a) 

r = |4jHi 2 /*«+W s + i^i 2 ^/*). ( 38b ) 

/ 3 = |^(^/'» + ^^) + i(l^| 2 -W 2 )/*. (38c) 



Substituting these expressions into the EL equations (|34ap and (j34bl) one obtains 

= ^if z -^f +^-0, (39) 
= ^ 2 f z z +iJif z +f 2 = 0. (40) 

In particular, the EL equations for tpi and ip2 (and the conjugates of these equations) involve f z z and 
f z (and their complex conjugates) but not f z z (or its complex conjugate). The tangential projection f z z 
remains undetermined at this level. At an algebraic level, this fact is related to the identity e z = ie z x n. 

The combination ipi e^ >1 + 4>2 £^, 2 determines f z z : 

ik e^i + V-2 e^ 2 = 2f z z + ViT 1 + j 2 f 2 = . (41) 
Using expressions ([33]) for T 1 and T 2 in this equation and solving for f z z one finds 

r ' z = ~\ {^Wi + ^ 2 ^h) + \ d ' z ^ 2 + ^ ~ iVIm ^^ ■ 

Now, by substituting into this last equation the expressions for A 1 and A 2 given in Q33p one finally obtains 
fz = ~\ (i>i^ + ^ 2 ^) + " \ ^2 2 d s c + ^ 2 d g c) . (42) 



7 Components of the complex conjugate of the spinor, and ip2, are varied independently of i/>i and ip2- 
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Similarly f z is determined by the combination ipi e^ 2 — ip2 '■ 

i>ie^-^ 2 e^ = \iP\ 2 f z + $ 1 T 2 -^ 2 T 1 = 0. (43) 
Substitution of expressions for T 1 and T 2 into this equation and solving for f z gives 

+ 7^2 (A 2 a 2 ^ 2 - A 1 ^ + V(A> 2 + AVi)) • 

Using once again expressions (1331) in place of A 1 and A 2 along with the the expressions (|C.2I) and (|C3[) 
for the derivatives of if)\ and tp2i the following simplification results 

f = ~ d i T777777T7 + 7779 ^TT ~ & 



2|V-| 2 SVJ + |i/>| 2 V ^ 



-'2 



|0 | 2 -(l^i| 2 - l^| 2 Rr + ^xihdgc-^dnc) . (44) 

The only ambiguity in the components f z z and / z of the stress tensor is the one inherited from the solution 
ofEqj 



As noted previously, the Euler-Lagrange equations for the spinor and the potential leave completely un- 
determined the off-diagonal component f z z of the tangential stress. At first, this appears to suggest that 
something is amiss. One must remember, however, that by representing the surface isothermally, one 
necessarily foregoes access to reparametrization invariance. This feature manifests itself in the tangential 
projections of the conservation law for the stress tensor. Whereas these equations would be satisfied iden- 
tically in any completely reparametrization invariant framework, in this one they provide the differential 
equations determining the missing component of the tangential stress. 

Taking the projections of the conservation law £x = 0, where £x is given by Eq. (|2"T]) . onto the tangent 
vectors provide the equation 

ex • e 2 = 8- z f z + H 4 d 2 (t^tA) + Af z + H 2 Vr = 0, (45) 

along with its complex conjugate expression. Thus far, it has not been necessary to specify explicitly the 
functional form of L. To solve Ea. (|4"5|) . we will suppose for simplicity that L depends only on \ip\ 2 and V, 
undifferentiated. These differential equations can then be solved for the missing component of the stress 
tensor. The most general solution is given by 

A = Tfj^^ + #i«r + \ {^ 2 d- z c + 4, 2 d- z c) + h(z) , (46) 
where h(z) is an arbitrary function. 

The projection onto the normal vector provides the "shape" equation 

ex ■ n = d z f z + d- z f - | (f z z + f%) - {Af% + Af z ) = . (47) 

We are now in a position to examine the different ambiguities which have arisen in our construction of the 
stress tensor. The first of these originates in the solution of the Euler-Lagrange equation for the potential. 

The contribution to the stress originating in the homogeneous solution AJj is given by 

|V| 4 f c 2 = d z c^ 2 \ z - V^e. + IVlViV^n) 

+ d- z c [4> 1 2 e z - <A 2 2 e 2 - - M 2 V>iV>2n) • (48) 
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Decomposing the complex- valued function c into its real and imaginary parts, c = c x + ic v , it is possible 
to express this contribution as the partial derivative of a space vector 

f| = id$(cyx.i +c^x 2 ) , (49) 

where the fact that the basis vectors of E 3 are constant has been used (holomorphic will do). With respect 
to the surface adopted basis these three vectors are given by 

xi - ^(^ 2 2 ~^ 1 2 )e z + ^ 2 2 -^ 1 2 )e- z + \^\ 2 (^ 2 + ^ 2 )n) 1 (50a) 

x 2 = ^(-(V'2 2 +V i i 2 )e 2 + (^ 2 2 +^ 1 2 )e 2 --|^| 2 (V'iV'2-^i^2)n), (50b) 

x 3 = ^ ^ 2 e 2 + ^2e 2 + ^| 2 (|V'i| 2 -|^2| 2 )nJ ■ (50c) 

Similarly, the contribution to the stress tensor originating in is given by 

f r 2 = id- z (rx 3 ) , (51) 

so that in full the contribution to the stress tensor arising from the homogeneous solution can be expressed 
as the partial derivative of a real- valued space vector V, 

fjf=tfcV, (52) 

where V is given by V = (cj,xi + c x ic 2 + rx^). Since V is real-valued, f z h has zero divergence: d z f z h -\-dzf z h = 
0. It is a null tensor which is automatically conserved and does not contribute to the shape equation (1471) . 
Therefore, as claimed above, it is legitimate to neglect this contribution to the stress tensor and retain 
only the part arising from the inhomogeneous solution A,;. 

Unlike the canonical stress tensor in the parametrized description of a surface, its counterpart in this 
framework is not unique. 

The ambiguity associated with the arbitrary function h(z) appearing in the solution of Eq. (|45l) given by 
Eq. (|46[) appears to be of a more serious nature. It is not simply a gauge artifact, contributing as it does 
to the shape equation a term 

2 

* W 4 

However, the normal projection of the Euler-Lagrange derivative ex • n associated with any reparametriza- 
tion invariant energy function had better form a scalar (density). Under a holomorphic transformation 
z — > w(z) the factor _4/|?/;| 4 transforms as 

A{z,z) w'(z) A(w,w) 

|V(z,r)| 4 w'(z)\^(w,iu)\' 1 ' 1 ' 

In order to form a scalar h(z) should transform as h(z) —¥ w' \z)/w' (z)h(w), which is not a holomorphic 
function. Consistency requires that h must vanish. 

In summary, for a functional which depends on \if)\ and V, the components of the stress tensor assume the 
simple form 

ft = l 54b > 

The Euclidean invariance of the energy and the identification of f 2 with the stress tensor within this frame- 
work is discussed in |Appcndix E| A more general dependence on the spinor and potential is considered 
in |Appcndix F 



4 (Ah(z) + Ah(z)) 
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3.1 Canham-Helfrich energy 

In soft matter applications, the surface energy will typically be a sum of bending energy, a term linear 
in K reflecting an asymmetry between the two sides of the surface, and an area term associated with a 
constraint or penalty on the total area [IB] , 

L = ~kL 2 +PLi+<jL 0) (55) 

where L n = l/2\^\ 4 K n . 
One finds that 

^ = (1-2)^^. and ^=2n\^K-\ (56) 
From equations (|54p one reads off the contribution of L n to the various components of the stress tensor: 

U\ = (57a) 

f n z - z = nAK n ~ 1 , (57b) 
/„* = -nd- z K n ~ x . (57c) 

It is simple to confirm that these expressions reproduce the well-known result |16) . 

in this particular parametrization. The corresponding contribution to the normal component of the Euler- 
Lagrange derivative (|47p is also easily shown to be given by 



£„ = -nAK™- 1 + K n ~ l (2nK G - (n - l)K 2 ) . (59) 

n=0: Area 

For surface area, corresponding to Lo = 1/2|?/>| 4 , only the diagonal tangential stress is non- vanishing: 

f z z = -\\M\ A = o, r = o. (60) 

The Euler-Lagrange derivative is given by £q — K; the critical points of area are minimal surfaces satisfying 
K = or V = 0. 

n=l: Integrated mean curvature 

For an energy proportional to the mean curvature, corresponding to L\ — l/2\ip\ i K , the components of 
the tangential stress tensor are 

A = -jH 4 ^, f'n = A, (61) 

and the normal stress vanishes, f z — 0. The Euler-Lagrange derivative £\ = 1Z; the critical points are 
developable with vanishing Gauss curvature. 

n—2: Canham-Helfrich bending energy Lb = Lij2 = 4V 2 
The components of the stress tensor are 

f', = 0, f' s =AK, f* = -e s K, (62) 
and the Euler-Lagrange derivative given by 

£ B = -AK + K (2K G - ^K 2 ^j . (63) 

The vanishing of f z z in this case is a manifestation of scale invariance. 



It was pointed out in the paragraph following Ea. ((25|) that if the constraints defining the stress tensor- 
relating the tangents to the embedding variables are not enforced, the Euler-Lagrange equations obtained 
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are incorrect. To underscore this point, consider the consequence of dropping this constraint so that 
f° = 0; the surface states are stress- free. In the case of pure bending, setting f a = in Eq. ([55|) for 
n = 2 implies that K = or K a b = g a t,K/2; the latter possibility implies a spherical geometry. The 
only stress-free solutions are thus minimal surfaces or spheres. Even though the variations of ip and V 
are consistent with the Dirac equation, and thus represent a surface, without the additional constraints, 
the corresponding Euler-Lagrange equations do not describe the critical points of the surface problem 
correctly. 

4 Discussion 

We have provided a variational framework, tailored to the representation of the surface geometry in terms 
of a spinor field interacting through a potential, to describe the equilibrium properties of this surface. The 
construction of a conserved stress tensor lies at the center of this framework. Our primary aim, of course, 
was not to provide yet another construction of the stress tensor. Even if this derivation does provide 
significant new insight into its relationship with the surface geometry, there are other derivations; rather 
it has been to examine how these variables can be reconciled with surface variational principles. One is 
now in a position to import the techniques of complex analysis to re-examine various inherently non-linear 
problems in soft matter where our toolbox has been found wanting. 

A problem that lends itself to be treated using the spinor representation of the surface geometry is a very 
old one that has been the focus of a revival of interest in recent years: how does a surface bend when 
its local geometry is constrained; a good example is provided by an unstretchable planar sheet of paper. 
Contrary to initial expectations, this is not a simple problem [29, 30J 131] ■ The general configuration will 
consist of piecewise developable surfaces (with vanishing Gaussian curvature) meeting along a set of ridges. 
Such surfaces assume a particularly simple form in the spinor representation on account of Eq. (IB. 41) . If the 
Gaussian curvature vanishes, \ip\ itself is harmonic. Unfortunately, developable surfaces do not generally 
occur as solutions of the unadorned Euler-Lagrange equations for bending: the correct Euler-Lagrange 
equations possess an additional term arising from the local constraint on the metric that needs to be 
imposed to maintain flatness under variation; they do not minimize bending energy unless one imposes 
this additional local constraint [32] . One also finds that the boundary conditions that are appropriate will 
reflect the constraint on the metric. A simple but non-trivial example is provided by a cone; it does not 
minimize bending energy unless the metric is constrained to be flat. This is why paper folds into cones but 
fluid membranes do not! This local constraint is very different in nature from the structural constraints 
appearing in the auxiliary framework; whereas the role of the latter set of constraints is to ensure that the 
spinors and the potential are consistent with some spatial geometry, they do not place any constraint on 
the geometry itself. It is remarkable therefore that the local constraint can be treated technically on an 
identical footing [35]. The extra term takes the form of a linear coupling of the corresponding tangential 
stress to extrinsic curvature. We are currently examining this problem within the spinor framework. 

There are several interesting directions for future work. Statistical field theory is an obvious one, where we 
have been largely limited to Gaussian approximations using traditional methods. There is, however, still 
some spadework to be done before one is in a position to do this with any confidence. It would be useful 
to first address a few questions of an elementary nature: how does one describe perturbation theory about 
some given equilibrium surface in this representation? Does the simple form of the bending energy have 
a counterpart in the expansion at second order? In this context, it is clear that even if one is interested 
in minimal surfaces described by the WE representation with V = 0, one needs to introduce a potential 
in the variation. Finally, it would appear that the coupling of electrons to the curved surface geometry in 
graphene calls naturally for a spinor description of the surface. 

Partial support from DGAPA PAPIIT grant IN114510 is acknowledged. 

Appendix A Rotation of the surface in E 3 

The Dirac equation is not SU (2) invariant. Nor are rotations of the surface realized in what would 
appear to be the "obvious" way by the action of SU(2) on the spinor ip using the local isomorphism 
50(3) = SU(2). Note, however, that the surface which corresponds to the spinor £ is the same as that 
obtained from ip under a rotation by an angle n around the X2 axis, i.e. under the substitution -0 — > £ the 
components of X change as X 1 —> ~X X ,X 2 — > X 2 ,X 3 —> —X 3 . More generally, the transformed spinor 
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ip = atp + where a and b are two complex numbers, also satisfies the Dirac equation so that it, as well, 
describes a surface X. If \a\ 2 + \b\ 2 = 1, ip is indirectly related to ip by an element U of SU(2) given by 

which can also be represented as J7 = e~ t ^ cr ' r , where f is a unit vector in E 3 and <x is the vector in E 3 with 
the Pauli matrices as its components. Because of the local isomorphism SO(3) = 577(2) this identifies 
the element R of SO(3) given by R = e 4,5J ' r , where J is a vector in E 3 with the infinitesimal generators 
of rotations as its components, so that R describes a anticlockwise rotation by an angle %9 around the 
direction f . Evidently ip ^ Uip, thus the manner in which U relates ip and ip is not through the usual 
action of SU(2), but with its associated rotation R that relates the embedding functions of both spinors, 
namely the embedding functions X describe a rotation of X, i.e. X = RX [12] , If \a\ 2 + \b\ 2 ^ 1, the 
transformation %p — ¥ aip + b£ describes a rotation accompanied by a scaling. 

Appendix B The Laplace-Beltrami operator and the Ricci scalar 

The Laplace-Beltrami operator A on the complex plane is given by 

A=p^ f . (B.l) 

The non-vanishing Christoffel symbols constructed with the metric g a b are 

Y\ z = g zs 8 z g z - z = d z \ng zz = 2 d z ln|^| 2 , (B.2) 

and its complex conjugate r z g2 . Likewise, the only non-vanishing components of the Riemann tensor in 
the conformal parametrization are 

R\- zz = ^A\n\^\ 2 = R s 2z - z . (B.3) 
The Ricci tensor (proportional to the metric) and the Ricci scalar are 

R zz = -^Aln\yj\ 2 , ft = -2Aln|V;| 2 . (B.L) 

Appendix C Identities for d z i/j\ and 

The two symmetric curvature scalars depend on ip only through the combinations \ip\ 2 and A. It is useful 
to possess identities for the missing partial derivatives, d z xp\ and d z ip 2 , in terms of these variables. 

Begin by differentiating \tp\ 2 with respect to z, using the fact that tp satisfies the Dirac equation, to obtain 

d z \ip\ 2 = ip 2 d z ip2 + i>id z ipi . (C.l) 
Now multiply across by -01 and add \p 2 A to both sides to get 

^A + ^dM 2 = M 2 0rfk, 

we thus obtain the identity 

d z ^ =^d z \n\rP\ 2 + -^rP 2 A. (C.2) 
Multiplying Eq. (|C.ip across by ip 2 and subtracting {pi A from both sides yields 

^ 2 dM 2 -xp~iA= \xP\ 2 d z ip 2 , 

giving the second identity 

d z 4>2=iP 2 d z \n\ip\ 2 --^rP 1 A. (C.3) 

m 

All derivatives of tp can now be expressed in the compact form, 

( d z ln\^\ 2 M~ 2 A\ ( V \ 
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Appendix D Gauss- Weingarten equations and the integrability 
conditions 

The Gauss- Weingarten equations, d a e a = r c ah e c — K a i,n and <9 a n = K a b e b , describe how the adopted 
frame {e a , n} changes as it is moved across the surface. In the generalized WE representation, by writing 
the adapted basis to the surface as E = {e z , e z , n}, these equations can be expressed in the compact form 

m 

<9 Z E = M E <9 Z E = NE, (D.l) 
where the linear transformations M and N are given by 

2<9 2 ln|V>| 2 -A \ / -jK\i>\ 4 

M=( -\K\^ , N= 0_ 2d z ln|V>| 2 -A I. (D.2) 



The integrability condition E zz = E zz leads to M z — N 2 + [M, N] = 0. From this condition and the linear 
independence of basis E one obtains the following relations 

**'»'*i 2 - w- v2j w J € K2 - (D - 3a) 

dr,A = \<l,\'d,V-Vd,\i,\ 2 = ^a.,K, (D.3b) 

which arc fulfilled by spinors and potential satisfying the Dirac-type Eq.JT]). So it is unnecessary to 
implement them in the effective functional when performing the variational principle. 
Equation (|D.3ap is the Gauss-Codazzi equation since it is equivalent to the equation 

U = K 2 - K ab K ab = K 2 - (K ZZ K ZZ + 2K zf K zS + K- z - z K zi ~). (D.4) 

Likewise, equation (|D.3b[) and its complex conjugate are the Codazzi-Mainardi equations for they are 
equivalent to the pair of equations 

V a if fcc = W b K ac or S7 Z K ZZ = W Z K ZZ and its c.c. . (D.5) 

Appendix E Conserved stress and torque 

The current associated with the variation of the embedding functions, denoted by Qx, is identified with 
the boundary term in the variation of the functional, H c , with respect to X, or 

Qx = -i fdzA dz (<9 2 (f 2 • SX) + d z {f z ■ 5X)) . (E.l) 

The current associated with the variations in the spinor field is given by 

= i J dz Adzd z (\ 2 5ipi), Q^=if dz A dzd^X'S^), (E.2) 

together with their complex conjugate expressions. 

The potential does not give rise to a current. Thus the complete variation of H c is given by 

SH c = i I dzAdze-SX + Q, (E.3) 



where the total current Q is given by Q = Qx + (Q^i + Qif>2 + c.c). When the Euler-Lagrange equations 
are satisfied, SH C = Q. 

Consider now a patch of surface bounded by a set of closed curves. To determine the change in the 
equilibrium energy under the translation of one of these curves (say T), consider a deformation SX that 
reduces to a translation S& on this curve and vanishes on the remaining boundaries. The spinor field 
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transforms trivially under translation. Making use of the Gauss theorem the change in the energy can be 
recast as 

SH C = -5a ■ F , (E.4) 
where the vector F is defined by the line integral 

¥ = J dsl z i z +c.c. . (E.5) 

with ds the line element along T, l z and its c.c. are the components of the covector associated with the 
vector of the Darboux frame adapted to T and which is normal to it but tangent to the surface. Eg. (IE. 4ft 
identifies the vector F as the force acting on the boundary curve T 24 . Furthermore, the quantity l z i z 
represents the local force acting on the line element ds, so that the conserved current associated with 
translational invariance P is correctly identified as the surface stress tensor. 



Rotational invariance is a little more involved due to the non-trivial transformation properties of the 
spinor field. Using the identity (1131) for the normal vector n in terms of the spinor field, it is easy to see 
that the variation in the spinor field induces a variation in n given by 

2 - - - - 

Sn = -— (Oi<5V>2 - ip2S^i)e z + (^1^2 - ^2Stpi)e z ) . (E.6) 

Making use of expressions (|3"3"|) for A and (|E.6I) for Sn in the expressions for the currents given above, it 
is easily seen that Q is given by 

Q = -i J dz A dzd z (f z • SX + c z ■ Sn) + c.c. , (E.7) 

where c z = ^f^e*. 

Considering now a infinitesimal rotation by a constant angle Soj. The embedding functions transform by 
<5X = Su> x X; the normal vector also rotates, Sn = 5u> x n. Thus the energy changes by 

SH c = —iSuj- J dz A dz d z m z + c.c. , (E-8) 

where m z = X x f 2 + l/2\ip\~ 2 SL/SVn x e 2 . Using an identical argument to the one used to identify the 
force on the boundary curve T, we identify the change in the energy associated with a rotation of this 
curve 

SH C = -Suj ■ M , (E.9) 

where 

M = / dsl z m z +c.c. . (E.10) 



We thus identify M as the total torque acting on this boundary. m z and its c.c. are identified as the 
components of the surface torque tensor [3^. Using the identity e z = ie z x n, the component of the torque 
tensor can be rewritten in the form 

j XT 



Appendix F Consistency of Euler-Lagrange equations 

At the end of section [5] it was noted that the Canham-Helfrich bending energy differs from the conformally 
invariant Willmore energy H3, defined by Eg . ((22)) by a topological energy proportional to the Gauss-Bonnet 
invariant. The corresponding Euler-Lagrange equations therefore coincide; indeed, the corresponding local 
stresses also coincide. It is instructive to demonstrate this explicitly using the functional form of H$ in 
terms of |*4|. 

In this case derivatives d z ipi and d z ip2 appear explicit through their Wronskian A and its c.c. The 
component f z g given by Eg. (|46l) is replaced by 

l '^ l '' +v ^m^)-^omi)- (E - 1) 
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For a functional which depends on derivatives of ip no higher than first, the relevant components of the 
stress tensor are 



r- - 5SF^ v (*ife-*i&)' ,E ' 2b) 



Note that for functionals involving higher order derivatives, the corresponding functional derivatives will 
contain extra terms arising from integration by parts. 

The required partial derivatives are 

Irw(^'- 2 w4 wkr-m A ^ (E - 3b) 

Making use of the Codazzi-Mainardi Eq. (|D.3b|) we have that the corresponding functional derivatives are 

^ = -^d,K, ^-=^d- z K, § = 0. (E.4) 

Substituting these identities into expressions (fl"2l and ((31]) for the components of the stress tensor repro- 
duces a stress identical to that for the Canham-Helfrich energy written in Eg. ([52]) . It is worth pointing 
out that the two stress tensors did not need to coincide: they could have differed by a null stress. 



References 

[1] Osserman R 1986 A Survey of minimal surfaces, second ed. (New York: Dover Publications Inc.) 

[2] Two very accessible introductions to minimal surfaces are provided by Gray A 2006 Modern Differen- 
tial Geometry of Curves and Surfaces with Mathematica (Boca Raton, FL: Chapman & Hall/CRC) 
and Fomenko A T and Tuzhilin A A 1991 Elements of the Geometry and Topology of Minimal Sur- 
faces in Three- Dimensional Space Translations of Mathematical Monographs Vol. 93, Amer. Math. 
Soc. 

[3] Schwarz A 1890 Gesammelte Mathematische Abhandlungen 1 Springer Verlag, Berlin 

[4] Tate MW, et al. 1991 Nonbilayer phases of membrane lipids Chemistry and Physics of Lipids 57 147 

[5] Kenmotsu K 1979 Weierstrass formula for surfaces of prescribed mean curvature Math. Ann. 245 
89; 2003 Surfaces with constant mean curvature Translations of Mathematical Monographs v. 221, 
Amer. Math. Soc. 

[6] Sullivan D 1989 The spinor representation of minimal surfaces in space Unpublished notes 

[7] Kusner R and Schmidt N 1993 The spinor representation of minimal surfaces in space GANG Preprint 
111.27 

[8] Konopelchenko B G 1996 Induced surfaces and their integrable dynamics Stud. Appl. Math. 96 9 

[9] Kusner R and Schmidt N 1996 The spinor representation of surfaces in space GANG Preprint IV. 18 

[10] Konopelchenko BG and Taimanov I A 1996 Constant mean curvature surfaces via integrable dynam- 
ical system. J. Phys. A 29 1261 

[11] Bobenko A I 2008 Exploring Surfaces through Methods from the Theory of Integrable Systems: 
The Bonnet Problem, in Surveys on Geometry and Integrable Systems, Advanced Studies in Pure 
Mathematics, Mathematical Society of Japan, Tokyo, Vol 51, p. 1; arXiv:math.DG/9909003 



16 



[12] Taimanov I A 2006 Two-dimensional Dirac operator and surface theory Russian Mathematical Surveys 
61 79; arXiv:2006 math.DG/0512543 

Seifert U 1997 Configurations of fluid membranes and vesicles Advances in Physics 46 13 

Canham P 1970 The minimum energy of bending as a possible explanation of the biconcave shape of 
the human red blood cell J. Theor. Biol. 26 61; 

Hclfrich W 1973 Elastic properties of lipid bilayers: theory and possible experiments Z. Naturforsch. 
C 28 693 

Guven J 2004 Membrane geometry with auxiliary variables and quadratic constraints J. Phys. A: 
Math, and Gen. 38 L313 

Capovilla R and Guven J 2002 Stresses in lipid membranes J. Phys. A: Math, and Gen. 35 6233; the 
construction was previously presented in a relativistic context in Arreaga G, Capovilla R and Guven 
J 2000 Noether currents for bosonic banes Annals of Phys 279, 126 

Steigmann D J 1999 Fluid Films with Curvature Elasticity Arch. Rational Mech. Anal. 150 127 

Lomholt M A and Miao L 2006 Descriptions of membrane mechanics from microscopic and effective 
two-dimensional perspectives J. Phys. A: Math, and Gen. 39 10323 

Evans E 1974 Bending Resistance and Chemically Induced Moments in Membrane Bilayers Biophys. 
J. 14 923 

Hoffman D A and Karcher H 1997 Complete embedded minimal surfaces of finite total curvature 
Geometry V. Encyclopaedia of Math. Sci. Vol 90 p. 5 (Springer- Verlag, Berlin, Heidelberg) 

Kusner R 1988 Global geometry of extremal surfaces in three-space. PhD thesis, (University of Cali- 
fornia, Berkeley); 1991 Bubbles, Conservation laws and balanced diagrams Geometric Analysis and 
Computer Graphics in Mathematical Sciences and Research Institute Publications, edited by P. Con- 
cus, R. Finn, andD. Hoffman Springer, New York, Vol. 17. 

Kim K S, Neu J and Oster G 1998 Curvature-mediated interactions between membrane proteins 
Biophys. J. 75 2274 

A. Agrawal and D.J. Steigmann 2009 Modeling protein-mediated morphology in biomembranes 
Biomechanics and Modeling in Mechanobiology 8, 371-79. 

Miiller M M, Deserno M and Guven J 2005 Geometry of surface mediated interactions Euro. Phys. 
Lett. 69 482; 2005 Interface-mediated interactions between particles: A geometrical approach Phys. 
Rev. E 72 061407 

Deserno M 2009 Membrane elasticity and mediated interactions in continuum theory, in Biomedical 
Frontiers: Nano structures, Models, and the Design of Life, Handbook of Modern Biophysics, edited 
by R. Fallcr, T. Jue, M. Longo, and S. Risbud Humana Press, New York, Vol. 2, p. 41-74 

Powers T R 2010 Dynamics of filaments and membranes in a viscous fluid Rev. Mod. Phys. 82 1607 

A. Agrawal and D.J. Steigmann 2009 Boundary-value problems in the theory of lipid membranes 
Continuum Mech. Thermodyn. 21 57 

Willmore T J 1982 Total Curvature in Riemannian Geometry (Chichester: Ellis Horwood) 

Ben Amar M and Pomcau Y 1997 Crumpled Paper Proc. R. Soc. Lond. A 453, 729 

Cerda E and Mahadevan L 1998 Conical surfaces and crescent singularities in crumpled sheets Phys. 
Rev. Lett. 80 2358; 2005 Confined Developable Elastic Surfaces: Cylinders, Cones and the Elastica 
Proc. R. Soc. Lond. A 461 671 

Witten T A 2007 Stress focusing in elastic sheets Rev. Mod. Phys. 79 643 

Guven J and Miiller M M 2008 How paper folds: bending with local constraints J. Phys. A: Math 
and Theor. 41 055203 

[33] Miiller M M, Deserno M and Guven J 2007 Balancing torques in membrane-mediated interactions: 
exacts results and numerical illustrations Phys. Rev. E 76 011921 



17 



